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Abstract 

A  mathematical  model  of  the  three-dimensional 
two-phase  reacting  flows  in  gas  turbine  combustors 
has  been  developed  which  takes  into  account  the 
mass,  momentum  and  energy  coupling  between  the 
phases.  The  fundamental  equations  of  motion  of  the 
droplets  are  solved  numerically  in  a  Lagrangian 
frame  of  reference  using  a  finite  difference  solu¬ 
tion  of  the  governing  equations  of  the  gas.  Well 
known  relations  are  used  to  model  the  heat  and  mass 
transfer  processes  and  the  initial  droplet  heat-up 
is  allowed  for.  The  entire  fuel  spray  is  construc¬ 
ted  using  a  finite  number  of  size  ranges  obeying  a 
two  parameter  droplet  size  distribution.  The  results 
are  found  to  be  in  close  agreement  with  experimental 
data.  An  important  feature  of  this  analytical  tech¬ 
nique  is  that  it  permits  the  rational  selection  or 
speci f ication  of  fuel  nozzle  design. 

Introduction 

The  advent  of  more  powerful  digital  computers 
has  provided  the  means  whereby  mathematical  model¬ 
ling  can  be  applied  to  gas  turbine  combustor  pro¬ 
blems  to  facilitate  greatly  the  arduous  initial 
stages  of  combustor  design.  This  is  now  of  particu¬ 
lar  interest  in  view  of  the  current  demands  which 
combustion  designers  are  required  to  meet  in  parti¬ 
cular  efficiency  of  combustion  at  a  wide  range  of 
operating  conditions  and  the  strict  control  of  poll¬ 
utant  emissions.  The  latter  have  become  increasingly 
stringent  in  recent  years  for  economic  and  political 
reasons.  The  present  trend  is  away  from  the  tradi¬ 
tional  cut-and-trv  methods,  which  are  expensive  and 
do  not  necessarily  produce  the  optimum  answer, to¬ 
wards  fundamental  modelling  of  the  physical  and 
chemical  processes  occurring  within  the  combustor. 

In  this  study,  the  problem  of  three-dimensional 
two-phase  combustion  is  approached  with  the  aim  of 
producing  an  algorithm  based  on  fundamental  princi¬ 
ples  which  can  correlate  all  the  details  of  combust¬ 
ion  occurring  within  a  gas  turbine  combustor  can. 

It  has  been  shown  by  Mellor1  that  a  key  factor  in 
gas  turbine  combustion  is  the  effective  fuel  distri¬ 
bution  which  is  determined  by  the  fuel  preparation 
technique.  It  has  also  been  indicated  by  Faeth  and 
co-workers^,  that  the  behaviour  of  a  liquid  spray 
flame  approaches  that  of  a  gas  diffusion  flame  pro¬ 
vided  a  sufficiently  fine  spray  is  considered.  It 
is  therefore  essential  that  the  model  incorporates 
an  accurate  representation  of  the  droplet  size  dis¬ 
tributions  encountered  in  gas  turbine  combustors, 
and  solves  the  relevant  equations  for  trajectory 
and  evaporation  taking  into  account  the  heat,  mass 
and  momentum  interaction  between  the  phases. 

In  liquid-fuelled  gas  turbine  combustors  the 
liquid  disintegrates  and  atomizes  soon  after  leaving 
the  injector  forming  a  spectrum  of  droplet  sizes. 

The  droplets  evaporate  along  the  trajectories  that 


they  follow  leaving  behind  pockets  of  fuel  vapour 
which  later  mix  with  the  oxidiser  and  burn.  Al¬ 
though  a  large  body  of  literature  exists  on  the 
combustion  of  single  drops  of  fuel11  it  has  been 
recently  indicated  by  Chigier1*  that,  in  the  light 
of  recent  experimental  research,  the  validity  of 
this  information  for  spray  combustion  is  question¬ 
able.  There  is  evidence  that  the  spray  acts  as  a 
permeable  surface  of  fuel  vapour  the  properties  of 
which  depend  on  the  interaction  between  the  gas 
flow  and  the  individual  droplets. 

Previous  studies  by  Swithenbank  and  co-workers 
on  modelling  combustors  involved  the  adaptation  of 
a  finite-domain  procedure  for  the  computation  of 
gas  flows  developed  at  Imperial  College*’-^.  In  order 
to  model  the  liquid  spray  several  size  ranges  were 
each  treated  as  a  chemical  species  and  defined  in 
terms  of  an  elliptic  equation  in  the  solution  pro¬ 
cedure.  The  trajectories  of  the  individual  droplets 
were  analysed  in  two  dimensions  by  stipulating  that 
flow  in  the  vicinity  of  the  atomiser  is  axially 
symmetric  in  a  supplementary  study10. 

In  order  to  more  realistically  model  the  com¬ 
plex  interactions  between  the  two  phases  a  new 
strategy  has  now  been  adopted  in  which  the  gas  phase 
is  calculated  in  an  Eulerian  frame  of  reference 
whilst  a  separate  calculation  is  performed  for  the 
liquid  phase  in  a  Lagrangian  frame.  The  latter  is 
achieved  by  applying  the  Particle-Source-In-Cel 1 
method11  to  the  three-dimensional  geometry  of  the 
present  problem.  This  approach  has  several  advant¬ 
ages  over  the  previous  method  and  it  is  believed 
that  ultimately  this  will  enable  the  designer  to 
evaluate  a  particular  combustor  design  without  the 
need  to  build  a  prototype,  and  to  give  him  more  de¬ 
tailed  information  on  performance  than  is  at  present 
available . 

In  particular,  this  technique  reveals  details 
of  the  three  dimensional  spray  trajectories  and  the 
consequent  evaporated  fuel  distribution,  thus  per¬ 
mitting  the  rational  selection  of  fuel  nozzle  char¬ 
acteristics.  It  can  therefore  fill  an  important  gap 
in  current  combustor  design  procedure. 

Equations  to  be  solved 

Since  most  gas  turbine  combustors  exhibit 

a  degree  of  cylindrical  cyclic  symmetry  in  their  de¬ 
sign,  it  is  convenient  to  present  the  governing 
equations  in  cylindrical  co-ordinates. 

The  mean  motion  of  the  gas  is  determined  by 
the  well-established  differential  equations  of  the 
conservation  of  mass  and  momentum. 

Mass  conservation 
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6-momentum 
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In  these  equations,  u,  v  and  w  are  the  velocity 
components  in  the  x,r  and  6  directions;  p  is  the 
pressure,  p  and  y  are  the  density  and  the  viscosity. 
It  can  be  seen  that  the  mass  conservation  equation 
contains  the  source  term  m"  which  is  due  to  the 
mass  addition  from  the  spray  to  the  gas  by  evapora¬ 
tion,  and  the  momentum  interaction  between  the  two 
phases  is  contained  in  the  terms,  Fx,  Fr  and  Fq.  The 
usual  Reynolds  stresses  have  been  replaced  in  the 
above  equations  by  the  effective  viscosity  hypothe¬ 
sis,  where  y  represents  the  effective  viscosity 
(i.e.  the  molecular  viscosity  augmented  by  its 
’turbulent'  counterpart). 

The  distribution  of  the  turbulent  viscosity  is 
provided  by  the  k-E  turbulence  model  used.  It  en¬ 
tails  the  solution  of  two  transport  equations  for 
turbulence  characteristics,  namely  the  kinetic 
energy  of  the  fluctuating  motion  and  its  dissipa¬ 
tion  rate.  The  differential  transport  equations 
for  k  and  e  are 
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The  turbulent  viscosity  is  related  to  k  and  c  by 
dimensional  arguments  in  the  following  way: 


(8) 
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which  are  assigned  the  following  valves: 
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The  additional  conservation  equations  required 
for  a  chemically  reacting  system  are  those  for 
enthalpy  and  chemical  species.  The  equation  for 
stagnation  enthalpy  h  can  be  written  as 
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where  S^includes  radiation  sources,  shear  work  and 
spray-gas  interaction  terms,  and  is  the  turbulent 
Prandtl  number  for  enthalpy.  The  definition  of 
stagnation  enthalpy  is 

h  -  C  T  +  m,  H  *  i  (u2  +  v2  +  w2)  . (10) 

p  fu  2 

which  is  used  for  the  recovery  of  temperature.  In 
this  equation  mju  is  the  mass  fraction  of  unburned 
fuel  in  the  mixture,  H  is  its  heat  of  reaction, 

Cp  is  the  specific  heat  of  the  mixture  at  constant 
pressure  and  T  is  temperature. 


The  conservation  equation  for  a 
soecies  j,  takes  the  following  form 
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where,  mj  is  the  mass  fraction  of  species  k.  Rj  is 
the  mass  rate  of  creation  or  depletion  by  chemical 
reaction  and  S,  is  the  rate  of  creation  by  evapora¬ 
tion  from  the  liquid  phase.  If  it  is  assumed  that 
the  combustion  reaction  is  a  single  step  kinetically 
influenced  reaction  between  a  pair  of  reactants  to 
form  a  single  product,  it  is  necessary  to  solve  one 
such  equation  for  each  of  m£u  and  m0J(.  The  influ¬ 
ence  of  turbulence  on  the  rate  of  reaction  is  taken 
into  account  by  employing  the  eddy  break  up  model 
of  Spalding1^.  The  reaction  rate  for  mfu  is  in  this 
case  taken  to  be  the  smaller  of  two  expressions 
given  by  the  familiar  Arrhenius  formulation  and  the 
eddy  break  up  model.  The  latter  is  conveniently 
described  as 


fu.EBb 


-  CRPk 


1  3m,  2 

1  [(— ) 
3x  ’ 


3m,  2 

(-4^3 

or 


3m 


fu 


>  } 


Where,  Cr  is  a  constant  of  order  6. 


(12) 
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The  density  distribution  in  the  flow  is  pro¬ 
vided  by  the  equation  of  state; 
m. 

p  -  P/  RT  Z  7^  . (13) 

j  Mj 

Here  M;  is  the  molecular  weight  of  species  j  and 
R  is  tne  gas  constant. 


13,14 


The  effects  of  radiation  are  accounted  for  in 
the  present  study  by  reference  to  the  fluxes  in 
the  positive  and  negative  x,  r  and  0  directions 
The  six  first  order  differential  equations  for 
these  fluxes  can  be  transformed  into  a  set  of 
three  second  order  equations  for  the  composite 
fluxes  Rx,  Pr  and  R0,  which  are  expressed  as, 
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where,  a  is  the  flux  absorption,  coefficient,  s  is 
the  scattering  coefficient,  E=cjt4  is  the  black  body 
emissive  power  at  the  fluid  temperature  and  0  is 
the  Stef an- Boltzmann  constant.  The  contribution  of 
radiation  to  the  source  term  of  the  enthalpy  equa¬ 
tion  is  calculated  from 
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and  the  equations  of  trajectory  are 
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-r~  =  u  ,  —  =  v  and  r-r-  =  w 
dt  p  dt  p  dt  p 


,  (22).  (23)  .  (2-1 


where,  u  is  the  molecular  viscosity  of  the  gas,  4; 
is  the  liquid  density,  Cp  is  the  drag  coefficient 
and  iip,  Vp  and  Wp  are  the  components  of  the  gas 
velocity  in  the  x,  r  and  6  directions  respectively. 
Re  stands  here  for  the  relative  Reynolds  numbc-r, 
which  is  defined  as, 

Re  =  DgD]up  -  u|/pg  .  (25) 


The  Reynolds  number  dependence  of  the  draf  coeffi¬ 
cient  in  sprays  is  given  by  the  following 
relations 
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h ,  rad 

Those  source  terms  of  the  above  equations  (1), 
(2),  (3),  (4),  (9)  and  (11)  which  arise  from  the 
interaction  between  the  gas  and  the  droplets  are 
obtained  from  the  solution  of  the  equations  of  the 
spray.  Formulations  of  spray  motion  in  both 
Eulerian1^  and  Lagrangianl6  frames  have  been  re¬ 
ported.  The  latter  method  is  more  suitable  for 
flows  where  the  spray  is  initially  confined  to  a 
thin  sheet17  and  has  therefore  been  adopted  in  the 
present  study. 

The  initial  size  distribution  of  the  spray  is 
divided  into  an  adequate  number  of  discrete  inter¬ 
vals  each  represented  by  an  average  diameter.  If  it 
is  assumed  Chat  the  size  distribution  is  of  the 
Rosin-Ramroler  type,  the  mass  fraction  of  drops  of 
diameter  greater  than  D  is  given  by 

Mp  •  expL~(D/D)  ]  . (18) 

where,  D  is  the  size  constant  and  n  is  the  size  dis¬ 
tribution  parameter.  It  should  be  emphasized  that 
the  use  of  the  Rosin-Rammler  distribution  is  not 
necessary  for  the  calculation,  however  our  experi¬ 
ments  have  shown  that  the  characteristics  of  most 
nozzles  can  be  presented  in  this  compact  and  con¬ 
venient  form. 

When  all  external  effects  except  the  drag 
force  are  neglected1®*,  the  equations  of  the  motion 
of  each  of  the  droplets  representing  the  size  groups 
which  constitute  the  spray  can  be  written  as; 


For  evaporating  droplets  an  additional  equation 
which  gives  the  rate  of  change  of  diameter  with 
respect  to  time  is  also  needed.  This  can  be  expresses 
as : 

Q 
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where  the  second  bracketed  term  on  the  right  hand 
side  accounts  for  the  forced  convective  augmentation 
of  evaporation  and  Cg  is  the  vaporization  rate  con¬ 
stant  whose  value  depends  on  the  physical  properties 
of  the  surrounding  medium  as  well  as  the  fuel  itself 
A  widely  quoted  expression  for  Cj,,  derived  from  the 
quasi-steady  analysis  of  droplet  combustion  is  due 
to  Wise  and  Agoston2®1. 
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Where,  A  and  Cp  are  the  thermal  conductivity  and  the 
specific  heat  at  constant  pressure  of  the  surrounding 
gas,  L  is  the  latent  heat  of  evaporation,  and  Too  and 
Tl  are  the  gas  and  droplet  temperatures  respectively. 


To  allow  for  the  'heating-up'  time  of  the 
droplets  after  entry  into  the  gas  stream,  it  is 
assumed  that  evaporation  begins  when  the  drop 
temperature  reaches  the  boiling  point.  The  equation 
for  the  time  rate  of  change  of  drop  temperature  can 
be  expressed  as 

■  6A(2*0.:  Re^Pr3)  - — —  .  (29) 
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where  the  Nusselt  number  h^j  been  replaced  by  the 
familiar  Ranz  and  Marshall*-  correlation. 

The  location,  velocity,  size  and  temperature  of 
each  size  range  of  droplets  are  determined  as 


(19) 


The  location,  velocity,  size  and  temperature 
of  each  size  range  of  droplets  are  determined  as 
functions  of  time  from  a  set  of  the  above  ordinary 
differential  equations  and  auxiliary  algebraic  re¬ 
lations  subject  to  initial  conditions  at  t=0.  The 
droplet  source  terms  are  obtained  by  calculating 
what  is  lost  or  gained  by  the  spray  in  terms  of 
mass,  momentum  and  energy  as  the  droplets  enter  and 
leave  volume  elements  of  finite  dimensions  in  three 
dimensional  space.  It  is  easy  to  show,  for  example, 
that  the  mass  deposited  in  an  arbitrary  control 
volume  by  droplets  of  a  given  size  group  is  the 
product  of  the  nur'oer  flow  rate  and  the  difference 
between  droplet  mass  at  entry  and  at  exit.  The 
total  mass  addition  is  obtained  by  summation  over 
all  the  size  groups  whose  trajectories  pass  through 
the  control  volume  in  question 
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where  Q0is  the  mass  flow  rate  of  drops  of  diameter 
group  Dq  at  the  point  of  spray  sheet  break-up.  The 
enthalpy  and  momentum  interaction  terms  are  derived 
by  similar  reasoning. 

Solution  Procedure 


The  previous  governing  partial  differential 
equations  for  the  conservation  of  mass,  momentum, 
energy  and  chemical  species  for  the  gaseous  phase 
are  reduced  into  their  finite-domain  analogues  by 
integration  over  the  computational  cells  into 
which  the  combustor  is  divided-2.  All  the  dependent 
variables,  with  the  exception  of  velocity  compo¬ 
nents  are  calculated  and  stored  at  the  nodal  points 
which  these  cells  encompass.  The  velocity  compo¬ 
nents  lie  on  the  cell  boundaries. 


The  resulting  algebraic  equations  can  he  re¬ 
presented  in  the  following  common  form: 
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where  the  A's  are  the  coefficients  which  contain 
the  contributions  from  the  convective  and  diffusive 
fluxes,  Sy  and  Sp  are  the  components  of  the  linear¬ 
ized  source  term. 

In  the  near  wall  regions,  the  usual  wall 
functions  are  matched  with  the  finite-domain  equa¬ 
tions  to  preclude  fine  grid  calculations  in  this 
region.  The  usual  practice  is  to  cut  the  link  be¬ 
tween  the  boundary  and  near  wall  points  by  setting 
the  appropriate  coefficient  equal  to  zero,  and  to 
insert  the  wall  influence  by  way  of  the  linearized 
source  terms2^. 

The  set  of  simultaneous  algebraic  equations 
can  be  solved  by  a  semi-implicit  iterative  scheme 
which  starts  from  given  initial  conditions  for  all 
the  variables  and  converges  to  the  correct  solution 
on  the  completion  of  a  number  of  iterations,  pro¬ 
vided  that  the  spray-gas  interaction  terms  these 
equations  contain  are  supplied  from  the  solution 
of  the  ordinary  differential  equations  which  des¬ 
cribe  the  motion  of  the  droplets. 

Each  iteration  performs  the  following  steps: 


i)  The  u,v  and  w  momentum  equations  are  solved 
sequentially  with  guessed  pressures. 

ii)  Since  the  velocities  at  this  stage  do  not 
satisfy  the  mass  continuity  equation  locally,  a 
"Poisson-type"  equation  is  derived  from  the  con¬ 
tinuity  equation  and  the  three  linearized  momentum 
equations.  This  pressure-correction  equation  is  then 
solved  for  corrections  to  the  pressure  field  and 
consequent  adjustments  are  made  to  the  velocity  com¬ 
ponents. 

iii)  The  k  and  e  equations  are  solved  using  the 
most  recent  velocity  field  and  the  effective  vis¬ 
cosity  is  obtained. 

iv)  The  iteration  is  completed  upon  the  solution 
of  enthalpy,  species  conservation  and  radiation 
equations . 

The  sets  of  simultaneous  ordinary  differential 
equations  for  each  of  the  droplet  size  ranges  are 
integrated  numerically  by  the  Runge-Kutta  method  at 
suitable  intervals  within  the  above  described  iter¬ 
ative  solution  procedure  to  calculate  and  up-date 
the  distributions  of  the  spray-gas  interaction  terms 
The  gas  velocities  are  assumed,  for  this  purpose, 
to  remain  constant  within  their  respective  cel  Is. At 
the  end  of  each  Runge-Kutta  step,  the  instantaneous 
locations  of  the  droplets  are  checked  against  the 
boundaries  of  the  volume  elements  that  they  are  tra¬ 
versing  and  when  migration  into  a  neighbouring  cell 
occurs,  the  corresponding  values  of  the  gas  veloci¬ 
ties  are  substituted  into  the  equations  before  pro¬ 
ceeding  into  the  next  time  step.  This  process  con¬ 
tinues  until  the  droplet  disappears  by  evaporation 
or  impingement  on  the  wall. 

Resul ts 


A  simple  gas  turbine  combustor,  on  which  ex¬ 
perimental  data  is  available,  v  .s  used  to  validate 
the  mathematical  model  presented  in  the  preceding 
sect  ions . 


Fig.l.  (las  turbine  combustor 

The  geometry  of  the  combustor  can,  which  was 
designed  as  a  research  combustor,  is  shown  in  Fig.l. 
The  air  stream  through  the  swirler  enters  the 
chamber  with  a  finite  swirl  velocity.  Primary, 
secondary  and  dilution  air  streams  are  introduced 
via  three  sets  of  injection  holes,  each  consisting 
of  six  equally  spaced  circular  orifices,  giving  rise 
to  the  three-dimensional  nature  of  the  problem.  The 
presumption  of  cyclic  repetition  of  the  flow  enables 
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the  combustor  can  to  be  represented  by  a  single  60° 
sector,  which  is  divided  into  a  27  x  18  x  7  element 
grid  in  the  x,  r  and  8  directions  respectively.  This 
particular  combustor  does  not  have  any  wall  cooling 
slots  and  it  was  operated  in  a  large  plenum  to 
eliminate  any  axial  components  in  the  radial  jet 
flows.  The  incoming  air  flow  was  divided  so  that 
7.8Z  of  the  total  flow  was  admitted  through  the 
swirler,  and  25.5%,  29.9%  and  36.8%  through  the 
primary,  secondary  and  dilution  jets  respectively. 

The  fuel  spray,  which  is  of  the  hollow  cone 
type  with  an  included  angle  of  60°,  is  introduced 
from  an  axially  located  fuel  nozzle.  The  Rosin- 
Rammler  parameters  used  in  equation  (18)  to  repre¬ 
sent  a  typical  fuel  spray  were  obtained  experimen¬ 
tally^  using  a  Malvern  particle  sizing  instrument. 
From  the  results  typical  values  of  D  *  60p  and 
n»2.2  were  taken.  A  total  of  20  droplet  size  ranges 
and  18  angular  injection  locations  within  the  60° 
sector  of  the  combustor  can  under  consideration  are 
used  to  construct  the  complete  spray  cone.  The  de¬ 
tails  of  the  operating  conditions  for  the  combustor 
are  summarized  below 

Total  air  flow  rate 
Total  fuel  flow  rate 
Fuel/Air  ratio 
Air  temperature  at  inlet 
Swirl  number 
Initial  velocity  of 
droplets 

Air  velocity  through 
injection  holes 
Air  velocity  through 
swirler 

A  vector  plot  showing  both  the  direction  and 
the  magnitude  of  velocities  in  the  combustor  for 
the  cold,  droplet-free  flow  is  given  in  Fig. 2. 


=  2.125  x  10  "  kg/s 
=  6.25  x  10-4  kg/s 
=  2.94  x  10“2 
=  351  K 
=  0.8 

=  20  m/s 
=  151  m/s 
=  27.6  m/s 


Fig. 2.  Velocity  vectors  in  the  combustor 

corresponding  to  the  cold,  droplet-free  flow 

It  can  be  noted  that  air  leaving  the  swirler  de¬ 
flects  sharply  outwards  due  to  the  centrifugal 
field  generated  by  the  inlet  swirl  and  forms  a 
well-defined  vortex  in  the  vicinity  of  the  primary 
jets.  All  the  jets  penetrate  deep  into  the  flow 
forming  relatively  stagnant  zones  in  the  wall  re¬ 
gion  between  jet  inlets.  Downstream  of  the  dilution 
jets,  the  relatively  large  radial  velocities  due  to 
the  incoming  jets  decay  substantially  and  the 


profile  attains  a  uniform,  axial  character  with  a 
maximum  located  on  the  center  line. 


In  Fig. 3.  the  velocity  vectors  in  the  combustor 
for  the  hot  flow  are  presented.  As  can  be  observed 
most  velocities  are  far  in  excess  of  those  encount¬ 
ered  in  the  cold  case.  Hotter  regions  of  the  flow 
domain  exhibit  high  velocity  values  associated  with 
high  temperatures  and  corresponding  low  values  of 
density.  The  recirculation  vortex  is  far  stronger 
in  the  hot  case  than  in  the  cold  one.  It  must  be 
mentioned  however  that  in  this  region  of  the  flow, 
the  spray  acts  as  a  constant  source  of  mass  and 
momentum,  enhancing  greatly  the  energy  content  of 
back  flow.  The  presence  of  combustion  has  also 


Fig. 3.  Velocity  vectors  in  the  hot  reacting  flow 


helped  the  velocity  fi.eld  between  the  primary  and 
secondary  jets  to  gain  more  uniform  and  axial  des¬ 
cription.  The  exit  velocity  profiles  tend  to  be 
more  uniform  and  the  secondary  and  dilution  jets  do 
not  penetrate  as  deep  into  the  flow  as  in  the  cold 
case . 


The  distribution  of  the  mass  fraction  of  fuel 
evaporated  but  unbumt  is  shown  in  Fig. 4.  and  as 
expected  shows  that  the  majority  is  found  very 
close  to  the  fuel  spray  with  the  effect  of  recir¬ 
culation  accounting  for  the  fuel-rich  region  near 
the  axis. 


Fig. 4.  Distribution  of  mass  fraction  of  unbumt  fuel 
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Knowledge  of  the  distribution  of  this  entity 
is  important  for  quenching  studies,  leading  for  ex¬ 
ample,  to  the  prediction  and  minimization  of  un- 
bumed  hydrocarbons  at  engine  idle  conditions. 


The  temperature  distribution  along  the  com¬ 
bustor  for  the  K=4  plane  is  displayed  in  Fig. 5.  The 
high  temperature  region  follows  the  hollow  spray 


Fig. 5.  The  distribution  of  isotherms  in  the 
combustor 


cone  with  a  cooler  core  caused  by  the  major  re¬ 
circulation  eddy.  High  temperatures  also  occur  in 
the  vicinity  of  the  combustor  walls  downstream  of 
each  air  injection  position  suggesting  that  com¬ 
bustion  may  be  occurring  in  these  regions.  The  low 
temperature  primary,  secondary  and  dilution  jets 
can  be  seen  entering  the  combustor  and  giving  rise 
to  cool  core  regions.  The  exit  temperature  profile 
which  exhibits  a  minimum  on  the  axis  and  shows  a 
maximum  in  the  neighbourhood  of  the  walls  is  con¬ 
sistent  with  the  measurements.  Such  a  profile  has 
been  obtained  previously  by  Turan^5  from  a  pre-mixed 
pre-vapourized  fuel  model  and  therefore  cannot  be 
solely  attributed  to  the  hollow  cone  fuel  spray. 

The  characteristics  of  this  particular  chamber 
geometry  are  to  blame  and  the  results  suggest  it 
should  be  redesigned  with  more,  smaller  air  entry 
ports  for  practical  use  in  an  engine. 


Fig. 6.  Distribution  of  the  turbulence  kinetic 
energy  (m2/s2)  in  the  combustor. 


Fig. 7.  Droplet  trajectories  for  a  60°  sector  in 
orthographic  projection 

fuel  nozzle.  For  comparison,  a  45°  included  angle 
spray  injected  into  the  hot  flow  is  depicted  in 
Fig. 8.  A  smaller  proportion  of  the  droplets  are 
seen  to  be  hitting  the  walls  of  the  combustion 
chamber.  Contrary  to  the  behaviour  of  the  previous 
spray,  the  droplets  are  seen  to  penetrate  into  the 
secondary  zone.  It  is  apparent  from  the  figure  that 
the  three-dimensional  flow  field  has  a  strong  in¬ 
fluence  on  the  trajectories  of  individual  droplets. 


Fig. 6.  shows  the  distribution  of  the  kinetic  energy 
of  turbulence.  The  most  noticeable  feature  here  is 
the  turbulent  core  near  the  axis  which  extends  for 
a  large  part  of  the  length  of  the  combustion 
chamber.  This  large  turbulent  zone  will  be  of  in¬ 
terest  for  any  subsequent  study  using  the  chemical 
reactor  network  approach. 

The  droplet  trajectories  computed  during  the 
execution  of  the  two  phase  flow  algorithm  can  be 
plotted  graphically  and  are  displayed  in  Fig. 7. 

Here,  the  spray  cone  for  a  single  60°  sector  is 
shown  by  means  of  orthographic  projection,  where  the 
small  circles  indicate  bum-out  locations.  It  can 
be  seen  that  the  droplets  deviate  from  the  nominal 
spray  cone  by  the  action  of  the  gas  flow  quite 
substantially.  Especially,  the  sudden  deflection 
those  droplets  which  encounter  the  primary  air  jet 
suffer  is  worth  noting.  A  large  number  of  droplets 
arc  found  to  be  hitting  the  combustor  wall,  although 
the  majority  of  the  fuel  is  evaporated  near  the 


Fig. 8.  Droplet  trajectories  for  45°  cone  angle. 

Cone! usions 

A  mixed  Lagrangian-Eulerian  two-phase  flow 
prediction  algorithm  has  been  formulated  and  applied 
to  the  case  of  combustion  in  a  gas  turbine  combustor 
The  model  predicts  not  only  the  flow  and  turbulence 
field,  temperature  and  species  concentration  dis¬ 
tributions,  but  also  the  trajectories  of  the  drop¬ 
lets  which  constitute  the  spray. 

Despite  the  complexity  of  the  looping  iterative 
nature  of  the  solution  procedure  employed  in  the 
present  study,  a  better  rate  and  degree  of  converge 
than  that  reported  by  Turan^S  was  achieved. 

The  computed  results  reproduced  all  the  pert¬ 
inent  features  of  the  hot  reacting  flow  in  gas 
turbines  and  were  found  to  be  in  accord  with  what 
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experimental  data  was  available. 

It  is  believed  that  a  stage  is  now  reached 
where  the  mathematical  modelling  of  gas  turbine 
combustors  presents  a  practical  alternative  to 
empirical  methods  during  the  early  steps  in  com¬ 
bustor  design.  It  must  be  mentioned  however,  that 
the  generality  of  the  computational  procedure  re¬ 
mains  and  it  is  applicable  to  any  other  gas-particle 
flow  problem. 

It  is  endeavoured  at  present  in  Sheffield  to 
introduce  a  more  fundamentally  based  turbulence 
model  into  the  present  algorithm  and  to  account  for 
the  turbulent  diffusion  of  particles  or  droplets  by 
means  of  a  stochastic  particle  tracking  method. 
Another  area  of  improvement  involves  the  refine¬ 
ment  of  the  droplet  evaporation  model  in  the  light 
of  the  recent  developments^. 

It  is  hoped  ultimately  to  be  able  to  improve 
the  performance  of  the  experimental  combustor  by 
means  of  modifications  based  on  the  predictions 
of  the  mathematical  model. 
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20  Abstract  A  Mathematical  model  of  the  three-dimensional  two-phase  reacting 
flows  in  gas  turbine  combustors  has  been  developed  which  takes  into  account  the 
mass,  momentum  and  energy  coupling  between  the  phases.  The  fundamental  equations 
of  motion  of  the  droplets  are  solved  numerically  in  a  Lagrangian  frame  of  refer¬ 
ence  using  a  finite  difference  solution  of  the  governing  equations  of  the  gas. 
Well  known  relations  are  used  to  model  the  heat  and  mass  transfer  processes  and 
the  initial  droplet  heat-up  is  allowed  for.  The  entire  fuel  spray  is  constructed 
using  a  finite  number  of  size  ranges  obeying  a  two  parameter  droplet  size  distri 
bution.  The  results  are  found  to  be  in  close  agreement  with  experimental  data.  A 
important  feature  of  this  analytical  technique  is  that  iL  permits  the  rational 
selection  or  specification  of  fuel  nozzle  design. 
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